The  Laplace  Transform  is  frequently  encountered  in  mathematics,  physics,  engineering  and 
other  fields.  However,  the  spectral  properties  of  the  Laplace  Transform  tend  to  complicate 
its  numerical  treatment;  therefore,  the  closely  related  “Truncated”  Laplace  Transforms  are 
often  used  in  applications. 

We  have  constructed  efficient  algorithms  for  the  evaluation  of  the  left  singular  functions 
and  singular  values  of  the  Truncated  Laplace  Transform.  Together  with  the  previously 
introduced  algorithms  for  the  evaluation  of  the  right  singular  functions,  these  algorithms 
provide  the  Singular  Value  Decomposition  of  the  Truncated  Laplace  Transform. 

The  resulting  algorithms  are  applicable  to  all  environments  likely  to  be  encountered  in 
applications,  including  the  evaluation  of  singular  functions  corresponding  to  extremely  small 
singular  values  (e.g.  1CV1000). 


On  the  Analytical  and  Numerical  Properties 
of  the  Truncated  Laplace  Transform  II. 


Roy  R.  Lederman',  Vladimir  Rokhlim, 
Technical  Report  YALEU/DCS/TR-1507 
May  29,  2015 


t  This  author’s  research  was  supported  in  part  by  the  ONR  grants  #N00014-11- 1-0718,  #N00014- 
10-1-0570  and  #N00014-14-l-0797  and  the  NSF  grant  #1309858. 

1  This  author’s  research  was  supported  in  part  by  the  ONR  grants  #N00014-10-l-0570  and 
#N00014-14-l-0797. 

Applied  Mathematics  Program,  Yale  University,  New  Haven  CT  06511 


Approved  for  public  release:  distribution  is  unlimited. 
Keywords:  Truncated  Laplace  Transform,  SVD. 


1  Introduction 


The  Laplace  Transform  £  is  a  linear  mapping  L2[0,  oo)  — >•  L2[ 0,  oo);  for  a  function  /  G  L2[0,  oo), 
it  is  defined  by  the  formula 

COO 

(£(/))  (u)  =  /  e-to/(i)df.  (1) 

Jo 

As  is  well-known,  C  has  a  continuous  spectrum,  and  C~l  is  not  continuous  (see,  for  example, 
0)-  These  and  related  properties  tend  to  complicate  the  numerical  treatment  of  C. 

In  addressing  these  problems,  we  find  it  useful  to  draw  an  analogy  between  the  numerical 
treatment  of  the  Laplace  Transform,  and  the  numerical  treatment  of  the  Fourier  Transform  J7; 
for  a  function  /  €  L1(R),  the  latter  is  defined  by  the  formula 

/OO 

(2) 

-OO 


where  oj  €  R. 

In  various  applications  in  mathematics  and  engineering,  it  is  useful  to  define  the  “Trun¬ 
cated”  Fourier  Transform  J-c  :  L2[— 1, 1]  — >•  L2[— 1,1];  for  a  given  c  >  0,  J-c  of  a  function 
f€L2[-l,l]  is  defined  by  the  formula 

(Tc(f))(co)  =  j\~^f(t)dt.  (3) 

The  operator  J-c  has  been  analyzed  extensively;  one  of  most  notable  observations,  made  by 
Slepian  and  Poliak  around  1960,  was  that  the  integral  operator  Tc  commutes  with  a  second 
order  differential  operator  (see  ng).  This  property  of  J-c  was  used  in  analytical  and  numerical 
investigations  of  the  eigendecomposition  of  this  operator;  for  example  in  mum  m  m  na 

151 . 

For  0  <  a  <  b  <  oo,  the  linear  mapping  Caj,  :  L2[a,  b]  — >•  L2[0,  oo)  defined  by  the  formula 

(CaM)  M  =  f  e~tuJf(t)dt,  (4) 

J  a 

will  be  referred  to  as  the  Truncated  Laplace  Transform  of  /;  obviously,  Ca,b  is  a  compact 
operator  (see,  for  example,  0)- 

The  Singular  Value  Decomposition  (SVD)  of  Ca ^  has  been  analyzed,  inter  alia,  in  [3]  and 
[5]|;  Bertero  and  Griinbaum  observed  that  each  of  the  symmetric  operators  (£a,b)*  °  £a,b  and 
£a,ft°(£a,&)*  commutes  with  a  differential  operator  (see  |5]).  Despite  3.  6l  3  2.. T j ,  13.3  T],  much 
more  is  known  about  the  numerical  and  analytical  properties  of  Tc  than  about  the  properties 
of  £«,&• 

We  have  constructed  algorithms  for  the  efficient  evaluation  of  the  SVD  of  Caj,-  In  [12]  we 
introduced  algorithms  for  the  efficient  evaluation  of  the  right  singular  functions  and  singular 
values  of  £a,fe-  In  this  paper,  we  introduce  an  efficient  algorithm  for  the  numerical  evaluation 


1 


of  the  left  singular  functions,  and  an  additional  algorithm  for  computing  the  singular  value. 
Additional  properties  of  the  Truncated  Laplace  transform  will  be  discussed  in  upcoming  papers. 

This  paper  is  organized  as  follows.  Section  fj2]  summarizes  the  various  standard  mathemat¬ 
ical  facts  and  simple  derivations  that  are  used  later  in  this  paper,  as  well  as  various  properties 
of  the  Truncated  Laplace  Transform.  Section  §[3]  contains  the  derivation  of  various  properties 
of  the  truncated  Laplace  transform  which  are  related  to  the  left  singular  functions;  the  prop¬ 
erties  are  used  in  the  algorithms  presented  in  the  following  section.  Section  §0  describes  the 
algorithms  for  the  evaluation  of  the  left  singular  functions  and  singular  values  of  the  Trun¬ 
cated  Laplace  Transform.  Section  ^contains  numerical  results  obtained  using  the  algorithms. 
Section  1}6]  contains  generalizations  and  conclusions. 


2  Preliminaries 


2.1  The  Laguerre  Functions 

In  this  section  we  summarize  some  of  the  properties  of  the  Laguerre  Polynomials,  Generalized 
Laguerre  polynomials,  and  Laguerre  Functions.  The  Laguerre  polynomial  Lk{x)  of  degree  k  >  0 
is  defined  by  the  formula 

=  (5) 

m= 0  v  ' 

The  Generalized  Laguerre  Polynomial  of  order  a  >  —  1  and  degree  k  >  0,  denoted  by  L^\x), 
is  defined  by  the  formula 


LG*)- EM)' 


m= 0 


k  +  a\  1 


k  —  mj  ml 


-x 


(6) 


The  Laguerre  Function  of  degree  k  >  0,  which  we  denote  by  &k,  is  defined  via  the  formula 

=  e~x/2Lk(x).  (7) 

The  following  well  known  properties  of  the  Laguerre  Polynomials  and  the  Generalized  La¬ 
guerre  Polynomials  can  be  found,  inter  aha,  in  [I]. 


Lak~\x)  =  LUx)-Lak_  i(x) 


(8) 


d_ 

dx 


Lk(x)  =  -4^ 


(9) 


xLk(x )  =  -(k  +  l)Lk+i(x)  +  (2 k  +  1  )Lk(x)  -  kLk_ i(x) 


(10) 
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(11) 


e~xtLk(x)dx  =  (t-  1  ft^-1 

As  is  well-known,  the  Laguerre  Polynomials  Lq(x),  L\(x), ...  form  an  orthonormal  basis  in 
the  Hilbert  space  induced  by  the  weighted  inner  product 

roo 

( f,g)=  e~xf(x)g(x)dx.  (12) 

■Jo 

It  follows  that  the  Laguerre  functions  <bo(x),  form  an  orthonormal  basis  in  the  stan¬ 

dard  L2[0,oo)  sense. 

By  (1TTD  and  (JT]), 

+  <13) 

2.2  The  Legendre  Polynomials 

In  this  subsection  we  summarize  some  of  the  properties  of  the  Normalized  Shifted  Legendre 
Polynomials,  derived  from  the  standard  Legendre  Polynomials.  We  define  the  Normalized 
Shifted  Legendre  Polynomial  of  degree  k  =  0, 1, ...,  which  we  will  be  denoting  by  P£,  by  the 
formula 


Pfc*(x)  =  Pk( 2x  -  l)V2fc  +  l;  (14) 

where  Pk  is  the  Legendre  Polynomial  of  degree  k;  the  standard  definition  of  the  Legendre 
Polynomials  can  be  found,  inter  alia,  in  [lj. 

As  is  well-known,  the  Legendre  Polynomials  form  an  orthogonal  basis  in  L2[— 1, 1],  but 
they  are  not  normalized;  a  simple  calculation  shows  that  the  Normalized  Shifted  Legendre 
Polynomials  Pf ,  P* , ...  form  an  orthonormal  basis  in  L2[0, 1]. 

By  (H2),  the  first  Normalized  Shifted  Legendre  Polynomial  is  a  constant, 

P*(x)  =  1.  (15) 


2.3  Singular  Value  Decomposition  of  integral  operators 

The  SYD  of  integral  operators  and  its  key  properties  are  summarized  in  the  following  theorem, 
which  can  be  found,  for  example,  in  [20] . 

Theorem  2.1.  Suppose  that  the  function  K  :  [c,d]  x  [a,b\  — >•  M  is  square  integrable,  and 
T  :  L2[a,  b]  — >  L2[c ,  d]  is  defined  by  the  formula 

0 T(f))(x)=  f  K(x,t)f(t)dt.  (16) 

J  a 
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Then,  there  exist  two  orthonormal  sequences  of  functions  uq,u\,...,  where  un  :  [a,  6]  — >•  R  and 
uo,ui,  where  vn  :  [c,  d]  — >•  R,  and  a  sequence  ao,cti,  •••  €  R,  where  a o  >  «i  >  ...  >  0,  such 
that 


un(t)f(t)dt  \  vn(x)  (17) 

for  any  f  E  L2[a ,  b\.  The  sequence  ao,  ai, ...  is  uniquely  determined  by  K. 

The  functions  uq,  u\,  ...  are  referred  to  as  the  right  singular  functions,  the  functions  vo,  v\, ... 
are  referred  to  as  the  left  singular  functions,  and  the  values  ao,a±,...  are  referred  to  as  the 
singular  values  of  the  operator  T.  Together,  the  right  singular  functions,  the  left  singular 
functions  and  the  singular  values  are  referred  to  as  the  SVD  of  the  operator  T. 

It  immediately  follows  from  Theorem  12.11  that 


Tilin')  —  OinVn , 

(18) 

T*(vn)  =  anun. 

(19) 

Observation  2.2.  The  right  singular  functions  uq,  u\,  ..  of  T  are  eigenfunctions  of  the  operator 
T*  o  T  and  the  left  singular  functions  vo,v\,...  are  eigenfunctions  of  the  operator  ToT*;  the 
singular  values  cio,  ot\, ...  of  T  are  the  square  roots  of  the  eigenvalues  of  T*  o  T  and  ToT*.  In 
other  words,  for  every  n  =  0, 1, ..., 

((T*  o  T)(un))  (t)  =  J  K(x,t)^J  K(x,t)un(t)dt 


^  dx  =  a^un(r)  (20) 


(T(f))  (x)  =  (  j 

_ n  W  a 


and 


((T  °T*)(vn))(0  =  /  /  K(x,t)vn(x)dx)dt  =  a2nvn(£) 


(21) 


Remark  2.3.  The  function  K  can  be  expressed  using  the  singular  functions  as  follows  (see 

ED]), 


oo 

K(x,t)  =  y 'vn(x)anun(t)  (22) 

n= 0 

and  it  can  be  approximated  by  truncation  of  small  singular  values  (also  see  [20]): 

p 

K(x,t )  ~  Y2vn(x)anun(t)  (23) 

n=0 
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2.4  The  Truncated  Laplace  Transform 

Definition  2.4.  For  any  pair  of  real  numbers  a,  b,  such  that  0  <  a  <  b  <  oo,  the  Truncated 
Laplace  Transform  Ca,b  is  the  linear  mapping  L2[a,b\  — >•  L2[0,oo),  defined  by  the  formula 

Mf))( u)=  /Vto/(t)dt,  (24) 

J  a 

Obviously,  the  adjoint  of  La^  is 

/*oo 

/  e-^g^du.  (25) 

Jo 

The  operators  Ca,b  and  (£a,b)*  are  compact,  the  range  of  (£a,i>)*  is  dense  in  L2[a,  6]  and  the 
range  of  Catb  is  dense  in  L2[0,  oo)  (see,  for  example,  0). 


2.5  The  SVD  of  the  Truncated  Laplace  Transform 

By  Theorem  12.11  there  exist  an  orthonormal  sequence  of  right  singular  functions  uq,u\,...  € 
L2[a,  b\,  an  orthonormal  sequence  of  left  singular  functions  vq,  v\,  ...  E  L2[0,  oo)  and  a  sequence 
of  real  numbers  ao ,  ai, ...  €  M  such  that 


(£a,b(f))  M  =  ^2an 

n= 0 


un(t)f(t)dt  vn(u j) 


(26) 


and  for  all  n  =  0,1, 


(27) 


and 


(£a,b)  (Vn)  —  &n'Ujri') 


(28) 


an  >  an+ 1  >  0.  (29) 

Remark  2.5.  The  multiplicity  of  the  singular  values  of  Ca  j,  is  one  (see  0);  in  other  words, 
for  all  n  =  0,1, ... 


Oin  >  Otn+l-  (30) 

Remark  2.6.  According  to  Observation  12.21  the  right  singular  functions  uo,ui, ...  of  £aib  are 
eigenfunctions  of  the  integral  operator  {£a,b)*  °  £ a,b  :  L2[o,  b]  — >•  L2[a,  6]  given  by  the  formula 

(((A,, b)*  o  Cafi)  (/))  (t)  =  f  e~ut  e~“sf{s)ds\  d^  =  £  fi^mds, 
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(31) 

and  the  corresponding  eigenvalues  of  {£a,b)*  °  C-a,b  are  a^,  a2, ...,  where  an  is  the  singular  value 
of  Catb  associated  with  the  right  singular  function  un.  In  other  words, 

(((£a,b)*  o  £a,b)  (un))  (t)  =  f  —^—un(s)ds  =  a2nun(t).  (32) 

J  a  t  ~r  S 

Similarly,  the  left  singular  functions  vn  of  Cab  are  eigenfunctions  of  the  integral  operator 
£-a,b  °  ( £a,b)*  ■  [0,  oo)  — >•  L2[ 0,  oo)  given  by  the  formula 


((£a,fcO  (Ca,b)*)  (g))  (uj)  = 


=  /  e 


—cut 


jHe-^OtOdp)  d t  =  JQ 


oo  -a(u+p)  _  -b(ui+p) 


0J  +  p 


-g(p)dp, 


(33) 


and  the  corresponding  eigenvalues  Ca,b  o  (£a,b)*  are  aft,  a2, ....  In  other  words, 

f  oo  -a(aj+p)  _  -&(w+p) 

((£a,b  °  {£a,b)*)  (Vn))  (w)  =  /  2 - ; - ^n(p)dp  =  «2Un(w). 

Jo 


IjO  +  p 


(34) 


2.6  The  differential  operators  Dt  and  Dw  associated  with  the  singular  func¬ 
tions  of  Ca,b 

In  this  subsection  we  summarize  several  properties  related  to  the  differential  operator  Dt, 
defined  by  the  formula 

(A(/))  (t)  =  ^  (V  -  a2)(b2  -  i2)^/^))  -  ~  a2) fit),  (35) 

where  /  €  C2[a,b\;  and  properties  related  to  the  differential  operator  Du,  defined  by  the 
formula 


d2 

duj2 


.2  ^ 

do;2' 


+  (a2  +  bz) 


d 

d  uj 


UJ 


do; 


+  (■ 


a2b2uj2  +  2a2)  /(w), 


(36) 


where  /  G  C4 [0,  oo)  n  L2[0,  oo).  For  a  derivation  of  these  properties,  see  [5]. 

Theorem  2.7.  The  differential  operator  Dt,  defined  in  \35 1).  commutes  with  the  integral  oper¬ 
ator  (Ca,b)*  °  C-a,b>  (specified  in  (EZJ))  in  L2[a,b\.  In  other  words, 

Dt  o  {(£afi)*  o  £aib)  =  i(£a,b)*  o  Cafi)  o  Dt  (37) 
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Theorem  2.8.  The  differential  operator  Du,  defined  in  IS6\).  commutes  with  the  integral  op¬ 
erator  Ca o  ,  (specified  in  $EW)  *n  [0,  oo).  In  other  words, 

£a,b  O  (Ca,b)*  O  Dcj  =  Du  O  £a  b  o  (Cajb)*.  (38) 


Theorem  2.9.  The  right  singular  functions  uq,  u\,  .. 
are  also  the  eigenfunctions  of  Dt. 


(defined  in  [27\))  of  Catb  (defined  in  [ 24\ )) 


Theorem  2.10.  The  left  singular  functions  vq,vi,  ...  (defined  in  { Wfy)  of  Ca,b  (defined  in  (24\)) 
are  also  the  eigenfunctions  of  Du. 


We  denote  the  eigenvalues  of  the  differential  operator  Dt  by  Xo,Xi,  ...,  and  the  eigenvalues 
of  the  differential  operator  Du  by  %q,x*,  ....  By  Theorem  12.91  the  singular  function  un  is  the 
solution  to  the  differential  equation 


^  f(t2  -  a2)(b2  -  t2)^un{t)  )  -  2(t 2  -  a2)un{t )  =  xnun(t), 


(39) 


and  by  Theorem  12.101  the  left  singular  function  vn  is  the  solution  to  the  differential  equation 


d2  /  2  d2  ,  A 
^  T  d^vt(-w)) 


+  ( a 2  +  h2) 


6_ 

da; 


_d_ 

da; 


vk(u) 


+  (—a2b2Lo2  +  2a2)vk(oj)  = 


(40) 


Remark  2.11.  The  singular  values  an  (defined  in  (1271)  1  of  the  integral  operator  Ca>b  are 
known  to  decay  exponentially  as  n  grows;  consequently,  the  direct  numerical  computation  of 
the  singular  functions  of  Ca ^  beyond  the  first  few  singular  functions  is  impossible. 

The  differential  operators  Dt  and  Du  are  advantageous  in  the  numerical  treatment  of  the 
singular  functions  un  and  vn  because  their  eigenvalues  are  well-separated. 


2.7  Properties  of  the  right  singular  functions  un 

In  this  section  we  presents  some  of  the  numerical  properties  of  the  right  singular  functions  un, 
a  more  detailed  discussion  of  these  properties  is  found  in  [12]. 

We  find  it  convenient  to  define  the  function  fin(x)  on  the  interval  [0, 1]  by  the  formula 

ifn{x)  =  Vb  —  a  un(a  +  (b  -  a)x).  (41) 

We  introduce  the  notation  hn  =  (/iq  ,  hf ,  ...)T  for  the  vector  of  coefficients  of  the  expansion 
of  the  function  ifn  in  the  basis  of  Normalized  Shifted  Legendre  Polynomials  (defined  in  (1141)1: 
where  the  element  hlf,  is  defined  by  the  formula 

K=  f  ifn(x)Pf{x)dx,  (42) 

Jo 
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so  that 


V’nO)  =  ^hlP*(x).  (43) 

fc=0 

Theorem  2.12.  The  vector  of  coefficients  hn  is  the  n  +  1-th  eigenvector  of  the  five-diagonal 
matrix  M ,  defined  by  the  formula 


]\,T  ,  _ _ (K  —  1)  K _ 

fc_2’fc  4v/2fc^3(2fc-l)x/2fc+T’ 


Mk~i,k  — 


k3(l+/3) 


_  (— 4— 6/3— 2fc/J(2+3/3)+fc2(7+12/3+2/32)  +  (2fc3+fc4)(7+16/3+8/32)) 

Mk,k  ~  2(2fc-l)(2fc+3)  ; 

Mk+1,k  = 


(44) 


Mi 


k-\-2,k 


5fe+iy2fc+3’ 
(fc+l)2(fc+2)2 
’  4v/2fc+T(2fc+3)\/2fc+5  ' 


where 


P  = 


2  a 


(45) 


b  —  a 

3  Analytical  apparatus 

3.1  Expansion  of  vn  in  the  basis  of  Laguerre  functions 

Suppose  that  g  is  a  smooth  function  in  L2[0,  oo).  Then,  g  can  be  expressed  in  the  basis  of 
Laguerre  functions  &k  (defined  in  Q);  let  r/  =  (r/o,  r/i ,  ...)T  be  the  a  vector  where 


COO 

Vk=  g(u)$k(u)  dw, 

Jo 

then  clearly  T]  is  the  vector  of  coefficients  in  the  expansion, 

OO 

5(w)  =  ^?7fc4>fc(u;). 


(46) 


(47) 


k= 0 


We  introduce  the  notation  r)n  =  (t/q,  rjf, ...)  for  the  vector  of  coefficients  of  the  expansion 
of  the  left  singular  function  vk  (defined  in  (1271)1  in  the  basis  of  Laguerre  functions;  where  the 
element  rff  is  defined  by  the  formula 


COO 

Vk=  Vn(uj)^k(uj)  dw, 

Jo 


(48) 


so  that 


Urdu;)  = 


)  = 


k= 0 


(49) 


3.2  A  matrix  representation  of  the  differential  operator  D w  in  the  basis  of 

The  purpose  of  this  subsection  is  to  express  the  differential  operator  Du  (defined  in  ifTGll )  in  the 
basis  of  Laguerre  Functions  &k  as  the  matrix  M  described  in  Leininald.il  Theorem  13.21  shows 
that  the  matrix  M  is  in  fact  a  five-diagonal  matrix;  and  Corollary  13.31  provides  the  relationship 
between  the  eigenvectors  of  M  and  the  left  singular  functions  vn  defined  in  (I27|). 

Lemma  3.1.  Let  g  be  a  smooth  function  with  the  expansion  rj  =  (go,r]i,  ...)T  specified  in  m- 


9(u)  = 

k= 0 

Suppose  that  if  =  D^^g),  with  the  expansion  c  =  (co,ci,  ...)T  such  that 

OO 

=  ^Cfc4>fc(  u). 

k= 0 

Then, 

c  =  Mr /, 

where  the  matrix  elements  Mjj~  of  M  are  defined  via  the  formula 
Mjk  =  j  (pw  ($*;))  (cj)dw, 

with  0  <  j ,k  <  oo. 

Proof.  By  the  linearity  of  the  differential  operator  Du  (defined  in  (l36l)L 


=  (pu(g)')  M  =  (Aj($fc))  H- 

k= 0 

Combining  m  and  ([MD, 

OO  oo 

^cfe4>fc( u)  =  y^?7fc  (pu($k)^  (w). 

k= 0  k= 0 

Now,  by  multiplying  both  sides  of  (1551)  by  and  integrating,  we  have 


'  OO  \ 

'Yh'Hk  (pi, ;($fc))  M  |  $j(w)du;. 


\k= 0 


By  linearity, 


c3  =  5Z  rlk  (y  j0  $i(w)  (^($fc))  (w)dw 


(50) 

(51) 

(52) 

(53) 


(54) 

(55) 

(56) 

(57) 
□ 
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Theorem  3.2.  For  any  k  >  0 


D, 


M  = 


+ 


(4a2  —  l)  (462  —  l)  (k  —  1  )k 
16 

k2  (16 o?b2  —  l) 

— ^ - 7 - L$k-1  (W) 


®k- 2 (w) 


k(k  +  l)  (-48 a2b2  -  4a2  -  462  -  3)  +  (-16a262  +  12a2  -  462  -  l)  , 

+  — - — - ^ (w) 


(k  +  l)2  (l6a262  —  l) 

+  - - ^ - -*fc+i  (w) 

(4a2  -  1)  (462  —  1)  (jfe  +  2)(jfe  +  l) 

16 


^k+2  (w)  , 


(58) 


where  is  the  Laguerre  function  defined  in  &■ 

In  other  words,  M  is  the  symmetric  five-diagonal  matrix  with  non-zero  entries  defined  by 
the  formulae 


Mk-2,k 

Mk-l,k 

Mk,k 

Mk+i,k 

Mk+2,k 


(4a2  —  1)  (462  —  1 )  (k—  1  )k 
16  ’ 

k 2 (16a2&2  —  1) 

4  ’ 

k(k  +  1)  (— 48a262  —  4a2  —  462  —  3)  +  (— 16a262  +  12a2  —  462  —  l) 

8  ’ 

(k  +  l)2  (l6a262  —  l) 

4  ’ 

(4a2  -  1)  (462  -  1  )(k  +  2 )(fc  +  1) 

16  ' 

(59) 


Proof.  By  the  definition  of  t)u  in  (l36l). 


(4,($*))  (x)  = 

=  +  +  +  (_a2&2a;2  +  2a2)$fc(a;) 

A  somewhat  tedious  derivation  from  ([60]).  using  identities  (j8j) ,  Q  and  (fT0|).  yields  (1581). 


(60) 

□ 
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Corollary  3.3.  Suppose  that  rf1  =  {r)Q,rf{, 
the  expansion  of  the  left  singular  function  vr 
ofM: 


..)  is  the  vector  of  coefficients  defined  in  \),  in 
defined  in  Wft);  then,  rf1  is  the  n+l-th  eigenvector 


Mr,n  =  X*nVn ,  (61) 

where  M  is  the  five- diagonal  matrix  m,  Xn  are  the  eigenvalues  of  the  differential  operator 
Du,  and  k  =  0, 1,  2, .... 

Proof.  By  (l40l).  vn  is  an  eigenfunction  of  Du,  with  the  eigenvalue  y*: 


( DujVk )  (u)  =  xtvk(u), 

so  that 

[Du  (u)  =  (DuVk)  M  = 

\  \k= o  /  /  k= o 

Therefore,  by  Lemma  13.11  we  obtain  (13.311. 


(62) 


(63) 


□ 


Remark  3.4.  If  the  operator  D u  is  represented  in  the  basis  associated  with  Hermite  polyno¬ 
mials,  rather  than  in  the  basis  of  Laguerre  functions,  the  resulting  matrix  M  in  Lemma  13.11 
and  Theorem  13. 21  is  seven- diagonal  and  of  certain  analytical  advantage.  This  approach  is  under 
investigation  and  will  be  reported  at  a  later  date. 


3.2.1  The  special  cases  £1/2, 7/2  an(i  £1/27,1/2 

An  inspection  of  formula  (1591)  immediately  indicates  that  there  are  two  special  cases  in  which 
the  second  off-diagonal  of  the  matrix  M  defined  in  (1591)  vanishes,  and  M  becomes  a  tridiagonal 
band  matrix. 

The  first  case  occurs  when  a  =  1/2.  The  substitution  of  a  =  1/2,  6  =  7/2  into  (1591)  yields 
Mk~i,k  =2  (72  -  1)  k2, 

Mk,k  =\  (-72  -  2  (72  +  1)  k2  -  2  (72  +  1)  k  +  1)  ,  (64) 

Mk+l,k  if)2  —  l)  ( k  +  l)2, 

and  zero  in  all  other  elements  of  M.  In  other  words,  the  differential  operator  associated  with 
the  Truncated  Laplace  Transform  £1/2, 7/2  yia  1|38[)  is  a  tridiagonal  band  matrix  in  the  basis  of 
Laguerre  Functions. 
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The  second  case  occurs  when  b  =  1/2.  The  substitution  ofa  =  b  =  1/2  into  (1591)  yields 


Mk-l,k 

Mk,k 

Mk+l,k 


(72  -  1)  k 2 

4^2  ’ 

(2  (72  +  l)  k2  +  2  (72  +  l)  k  +  72 


472 

(72-l)  (k  +  l)2 

472 


1) 


(65) 


and  zero  in  all  other  elements  of  M .  Again,  the  differential  operator  associated  with  the 
Truncated  Laplace  Transform  £1/27,1/2  via  (T38l)  is  a  tridiagonal  band  matrix  in  the  basis  of 
Laguerre  Functions. 


3.2.2  The  special  case:  C  !  ^ 

An  additional  special  case  occurs  when  \fab  =  1/2,  i.e.  in  the  case  of  the  Truncated  Laplace 
Transform  C  x  3/^.  Substituting  a  =  ,  b  =  ^  into  (|59l)  eliminates  the  first  off-diagonal  of 

Vf’T  v  ' 

M  and  yields  the  matrix 


M, 


(s) 


k—2,k 


M, 


W 

k,k 


M, 


(s) 


k-\-2,k 


(7  —  i)2(fc  —  i)jfc 
167 

(— 72  —  67  —  l)  k(k  +  1)  —  72 
87 

(7-  l)2(fc  +  l)(fc  +  2) 

167 


27  +  3 


(66) 


with  zero  in  all  other  elements  of  M^s\ 


3.3  The  “standard”  form  of  the  Truncated  Laplace  Transform 

As  shown  in  [3],  the  behavior  of  the  Truncated  Laplace  Transform  Ca  b  is  determined  by  the 
ratio 


7  =  b/a.  (67) 

The  following  lemma  summarizes  the  connections  between  the  SVD  of  the  Truncated  Laplace 
Transform  Caj,  and  the  SVD  of  the  Truncated  Laplace  Transform  C  x  ^ ;  we  will  refer  to  the 

latter  as  the  standard  form  of  the  Truncated  Laplace  Transform. 

Lemma  3.5.  Suppose  that  0  <  a  <  b  <  00,  and  that  un,  vn  and  an  are  the  n  +  1-th  right 
singular  function,  left  singular  function  and  singular  value  of  Ca>b,  respectively.  Suppose  that 
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un,  vn  and  an  are  the  n  +  l-th  right  singular  function,  left  singular  function  and  singidar  value 
of  the  standard  form  L yy,  respectively,  with  7  =  6/0. 


Then, 


un{t )  = 


2^70 


un(t/2^ya) 


vn(uj)  =  ^2^/7  a  vn(uja2y/i), 


and 


(68) 

(69) 


Qn  —  OLn» 


(70) 


Proof.  The  identities  are  readily  obtained  by  the  change  of  variables  t!  =  t/2a^y  in  (|32l).  the 
change  of  variables  u/  =  2a^pyoj  in  (13411 .  and  normalization  of  the  singular  functions.  □ 


We  denote  the  differential  operator  associated  with  C 
differential  operator  is  specified  by  the  formula 


1  Vi 


via  Theorem  12.81  by  Du ;  this 


_cP_ 

d  uj2 


,  72  +  l  d 

4y  dcu 


+ 


~hu2  +  ^) /Mj 


(71) 


which  is  a  special  case  of  formula  (1361). 
the  eigenvalues  of  Du  associated  with  C 
associated  with  Ca,b- 


The  following  lemma  specifies  the  relation  between 

!  ^  and  the  eigenvalues  of  Dw  (defined  in  ([551)1 

Vf  ’~2~ 


Lemma  3.6.  Suppose  thatx*n  the  n+l-th  eigenvalue  of  the  differential  operator  associ¬ 
ated  with  Ca,b,  and  that  x*n  the  n  +  1-th  eigenvalue  of  the  differential  operator  associated 
with  L _j__  ^/jt  where  7  =  6/ a.  Then, 


Xn  =  47«2X*n  =  4abx*  n. 


Proof.  The  relation  is  obtained  from  (13611  by  the  change  of  variables  u/  =  2a^yoj. 


(72) 

□ 


3.4  The  operator  C1 


In  this  discussion,  we  find  it  useful  to  introduce  the  transform  C7  :  L2 
we  define  C7  of  a  function  /  by  the  formula 


1  yT 
2^7’  2 


= >  L2  [—  1, 1]; 


(»)  =  7*/4  /  (t*/2/2)  . 


(73) 
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Remark  3.7.  A  simple  calculation  shows  that 


1 


(C-7(/))  (s)  ■  (C^g))  (s)ds 


4 

log  7 


f(t)g(t)dt. 


(74) 


The  following  lemma  provides  an  expression  for  the  function  C7 


the  Laguerre  function  defined  in  JJ]). 


,  with  <4>fc 


Lemma  3.8.  For  any  1  <  7  <  00, 


(s)  =  7S/4 


(75) 


Furthermore,  the  function  C1  I  C*  1  ^(3>fc)  I  is  even  or  odd,  depending  on  the  value  of  k: 

\  Vf ’“2“  / 

(c7  (^^(Sfc)))  (s)  =  (-1)*  (c-7  (^\_  ,£(**)))  (-«)■  (76) 

Proof.  The  identity  (1751)  is  obtained  by  substituting  (fT3l)  into  ([73]).  The  symmetry  property 
(1761)  follows  immediately  from  ([  7 5 1) .  □ 


Corollary  3.9.  For  1  <  7  <  00,  the  functions 


as  k  grows,  in  the  following  sense 


(*) 


decay  exponentially 


to 


7s/2  -  1 

7S/2  _)_  1 


1 

2 


2 

7s/2  +  1 


k 


(77) 


Furthermore,  if—  1  <  s  <  1,  f/ien, 


(s) 


2 

71/2  +  1 


(78) 


Proof.  By  d75l). 


7«/4 

7V2   1 

7s/2  +  1 

7s/2  +  1 

(79) 


so  that  for  all  1  <  7  <  00  and  s  €  R, 


7s/2  -  1 
7s/2  +  1 


(80) 


The  inequality  (1771)  follows  immediately  from  (]79l)  and  (l80l) .  Equation  (1751)  follows  immediately 
from  (1771) .  □ 


14 


3.4.1  The  function  C7{un) 


We  introduce  the  notation  Un  for  the  C7  of  the  n  +  1  right  singular  function  of  C 
defined  by  the  formula 


i  \/7  j  Un  is 

2^7’  2 


Un(s)  =  (C7(Un))  ( S ). 


(81) 


where  un  is  the  n  +  1-th  right  singular  function  of  the  operator  C  _j_  y=y  and  C7  is  dehned  in 

Vf’TT 

(ESD. 

The  following  lemma  summarizes  some  of  the  observations  found  in  [5j  and  reformulates 
them  using  the  notation  used  in  this  paper. 


Lemma  3.10.  The  function  Un  (defined  in  UH\))  is  an  eigenfunction  of  the  differential  operator 
Ds,  defined  by  the  formula 


( Ds  °  /)  (s)  =  (log(VT))  2  h2  +  1  ~  2y  cosh  (2s  log  (v^y)))  ^/(s) 
-  (il  cosh  (2s  log  (Vt))  +  \t 2  -  7~lj  /(»)■ 


(82) 


Since  the  differential  operator  Ds  is  symmetric  around  0,  the  function  Un  is  even  or  odd  in 
the  sense  that 


Un(s)  =  (-1  )nUn{-s). 


(83) 


3.5  Decay  of  the  coefficients 

Since  the  left  singular  function  vn  (defined  in  (1271)1  is  a  smooth  solutions  of  a  differential 
equation  (specified  in  (TIOl) ) ,  we  expect  the  coefficients  rff  (defined  in  (jj8j) )  in  the  expansion  of 
vn  to  decay  rapidly.  In  this  section  we  provide  an  estimate  for  the  actual  decay. 


Lemma  3.11.  Suppose  that  vn  is  the  n  +  1-th  left  singular  function  of  the  operator  C  1  75 
(defined  in\24\);  then 


i y/2 


—  Of n  r. - 

v  log  7 


1  - 


1  +  ^7 


(84) 


where  rjf  is  defined  in  |75|)  and  an  is  the  n  +  1-th  singular  value  of  C  1  ^ . 

2 


Proof.  By  (HH1)  and  fl27l). 


hk 


1  r  f  rb 

J  0  \Ja 


e  Utun(t)dt  4>fc(o;)du;. 


(85) 
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Changing  the  order  of  integration  and  using  (j25l) 


Vk  =  “n1 


i  Vi 
2^7’  2 


(86) 


A  simple  calculation  using  (1741)  and  (1811)  shows  that 


Vk  =  ar 


-l 


1  \/7 

2^’  2 


(s)ds. 


(87) 


By  the  Cauchy-Schwarz  inequality, 


— 


-1 


I  ^  (£4(s))2d.si 


/- 1 


cl 


•C  l  ypy 

2^7’  2 


($* 


(s)ds. 


(88) 


Now,  since  the  right  singular  function  un  is  normalized,  and  using  (1741)  and  (1811). 


\Vk\  <  a 


-l 


VlogT 


(89) 


Finally,  using  inequality  (f78|). 


\Vk  \  ^  « 


i V2 


%/logy 


1  - 


1  +  V7 


(90) 


□ 


Remark  3.12.  Lemma  13. 1 II  can  be  generalized  to  the  following  bound  for  the  coefficients  of 
the  expansion  of  the  left  singular  functions  of  the  Truncated  Laplace  Transform  Ca  &  in  the 
nonstandard  form;  the  proof  is  found  in  our  report  m-  The  coefficients  rff  in  the  expansion 
of  the  left  singular  function  vn  of  Ca y  are  bounded  in  the  sense  that 


\Vk  \  <  « 


-l 


./  2 

rySmax /2  _ 

V  logy 

sy^max /2  _|_  ^ 

(91) 


where 


f 

log  2a 

log  2b 

Smax  —  max  ( 

logy 

5 

logy 

(92) 


16 


3.6  Integral  of  the  function  un 

The  following  lemma  describes  the  relation  between  the  integral  un(t)dt,  where  un  is  the 
right  singular  function  (defined  in  (1271)1.  and  the  expansion  specified  in  (1431) .  The  expansion 
(H31)  is  associated  with  the  right  singular  function  un  via  (ITT!)  and  (H21)  and  it  is  studied  in  |12|. 

Lemma  3.13.  Suppose  that  un  be  the  n  +  1-th  right  singular  function  of  Ca,b,  then 

fb 

/  un(t)dt  =  Vb  —  ahfj.  (93) 

J  a 

where  Hq  is  the  first  coefficient  in  the  expansion  /[j£\)  of  the  function  ifn  defined  in  in  the 
basis  of  Legendre  Polynomials. 


Proof.  By 

K  =  [  'f’n{x)Po(x)dx,  (94) 

Jo 

Substituting  (fl5l)  into  (1M|)  yields 

K  =  [  ^n{x)dx.  (95) 

Jo 

Now,  substituting  fE34|)  into  (1M|)  with  the  change  of  variable  t  =  a  +  (6  —  a)x  yield  (1931).  □ 


3.7  A  relation  between  un,  vn ,  and  the  singular  value  an 

The  following  lemma  provides  the  relation  between  the  coefficient  hrf  (see  (I42|) ) ,  the  left  singular 
function  vn,  and  the  singular  value  an. 


Theorem  3.14.  Suppose  that  un,  vn  and  an  are  the  n  +  1-th  right  singular  function,  left 
singular  function  and  singular  value  of  La,b  and  suppose  that  hfi  is  the  first  coefficient  in  the 
expansion  defined  in  fffffy.  Then, 

fan 

an  =  Vb-a  — (96) 
vn{0) 

Proof.  Evaluating  both  side  of  CUD  at  0  yields 


Q!n’^n( 0)  —  (•£(!, b(^hi))  (0)  —  ®n^n(0)  —  /  Un(t)dt . 


Substituting  l93l  into  (|97T) . 

otnVn{d)  =  Vb-  a  hf. 


(97) 


(98) 


□ 
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4  Algorithms 

4.1  Computing  the  left  singular  function  vn 


In  this  section  we  introduce  an  algorithm  for  the  numerical  evaluation  of  vn(uj),  the  n  +  1-th 
left  singular  function  (defined  in  (l27l)l  of  Ca^  (the  operator  defined  in  (|24l)l. 

By  Theorem  12.101  the  function  vn(uj)  is  an  eigenfunction  of  the  differential  operator  D u 
defined  in  (1361):  by  Corollary  13.31  the  expansion  r]n  =  (rff,  r/”,  ...)T  (defined  in  fl49ll)  of  vn  in  the 
basis  of  Laguerre  functions  is  an  eigenvector  of  the  matrix  M  defined  in  ([59D. 

Based  on  Lemma  [33]  and  Lemma  13.61  it  is  sufficient  to  compute  vn(u),  the  n  +  1-th  left 

singular  function  of  the  Truncated  Laplace  Transform  in  the  standard  form  C  1  ^ ,  where 

Vf  ’“2“ 

7  =  b/a;  vn(t)  is  computed  from  vn(t)  using  (ISUll.  In  this  special  case  of  the  Truncated  Laplace 
Transform,  the  expansion  r\n  =  (t/q  ,  77”,  ...)T  of  vn(uj)  is  an  eigenvector  of  the  matrix  M ^ 
specified  in  ([Mil. 

Since  the  matrix  M ^  is  five-diagonal,  and  is  non-zero  only  in  even  rows  of  even  columns, 
and  on  odd  rows  of  odd  columns,  an  eigenvector  of  M^s'  must  vanish  at  either  the  even  or  odd 
positions. 

Therefore,  for  any  even  n  =  2m  where  m  >  0,  equation  (1491)  is  reformulated  as 


V2m(w)  =  y>rW’m*2fcM,  (99) 

k= 0 

where  r)even,m  =  (r/Qtera,m,  ,  ...)t  is  simply  the  even  numbered  elements  of  the  vector  r/2m; 

for  all  integer  k  >  0: 


2m  even.m 

V2k  =Vk 


(100) 


In  other  words,  r]even’m  is  the  m  +  l-th  eigenvector  of  the  tridiagonal  matrix  M^s,even\  obtained 
by  removing  the  odd  numbered  rows  and  columns  of  the  non-zero  elements  of  M^s'even> 

are  specihed  by  the  formula 

ny(s,even)  __  (7  —  l)2(2fc  —  1  )k 
1V1k-l,k  ~  g7 

^s,even)  _2  (-y2  -  67  -  l)  fc(2 fc  +  1)  -  y2  -  27  +  3  ^101^ 

k,k  g7 

ny(s,even)  _  (7  ~  l)2 (2k  +  l)(/c  +  1) 

1Uk+l,k  ~  g7 

Similarly,  for  any  odd  n  =  2m  +  1  where  m  >  0,  equation  ©  is  reformulated  as 


V2m+l(oj)  =  ^rfkdd'm<S>2k+l(u), 
k= 0 


(102) 
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where  r/odd,m  =  (r]Qdd’m ,  r]°dd,m ,  ...)T  is  simply  the  odd  numbered  elements  of  the  vector  r/2m+1; 
for  all  integer  k  >  0: 


r>2k+ 1  —  % 


2m+l  _  odd,m 

2k+l  ~  'Ik 


(103) 


In  other  words,  rjodd,rn  is  the  m  +  1-th  eigenvector  of  the  tridiagonal  matrix  M^s,odd\  obtained 
by  removing  the  even  numbered  rows  and  columns  of  M the  non-zero  elements  of  M^s,odd^ 
are  specified  by  the  formula 


(7  -  l)2k(2k  +  1) 


87 

2  (— 72  —  67  —  l)  (2 k  +  1  ){k  +  1)  —  72  —  27  +  3 


(104) 


87 


(7-  l)2(fc  +  l)(2£:  +  3) 


87 


Therefore,  the  algorithm  for  computing  the  left  singular  function  vn  of  the  Truncated  Laplace 
Transform  Ca,bi  where  n  =  2m  is  an  even  number: 

Step  1:  Compute  ifven’n ;  the  n+  1-th  eigenvector  of  the  matrix  M(s,even\  defined  in  (jlOljh 
Step  2:  Compute  the  function  vn  from  rjeven,n ,  using  the  expansion  specified  in  (1991). 

Step  3:  Obtain  vn  from  vn  using  (mi) . 

Similarly,  the  algorithm  for  computing  the  left  singular  function  vn  of  the  Truncated  Laplace 
Transform  where  n  =  2m  +  1  is  an  odd  number: 

Step  1:  Compute  ri°dd,n ,  the  n  +  1-th  eigenvector  of  the  matrix  M^s,odd\  defined  in  ()1041) . 
Step  2:  Compute  the  function  vn  from  7] odd’n ;  using  the  expansion  specified  in  (I102D. 

Step  3:  Obtain  vn  from  vn  using  (PTT1) . 

Remark  4.1.  For  computations  to  precision  e,  the  vector  if1  =  (^q  ,  77”,  ...)T  is  truncated  at 
K,  such  that  |?y^|  <C  e  for  all  k  >  K.  By  lemma  13.111  the  coefficients  rf£  decay  rapidly  as  k 
grows;  since  the  vectors  7 jeven>n  anc[  rjodd,n^  COntain  only  the  even  and  odd  elements  of 

77^),  respectively,  these  vectors  are  truncated  at  approximately  K/2  elements. 

The  actual  position  of  the  last  significant  coefficient  of  rjeven’n  and  r)odd,n,  larger  in  magnitude 
than  e,  is  given  in  Figure  [6]  and  Table  [3]  in  Section  for  several  combinations  of  7  and  n. 

4.2  Computing  the  singular  value  an 

Theorem  13.  Ill  provides  formula  (1961)  for  computing  the  singular  value  an  using  the  function  vn 
and  the  first  coefficient  of  the  vector  hn  (the  vector  associated  with  the  right  singular  function 
un  via  dUD  and  fl32D).  The  function  vn  is  computed  using  the  algorithm  in  Section  114.11  and 
the  vector  hn  is  the  n  +  1-th  eigenvector  of  the  matrix  M  (defined  in  (1441).  see  |12|) . 
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Remark  4.2.  The  numerical  difficulty  in  this  algorithm  is  that  the  first  element  Kq  of  the 
vector  hn  (see  (@2]))  is  of  the  order  of  an,  so  that  is  decays  exponentially  as  n  grows.  It  has 
been  shown  in  [jLj]  that  in  some  band  matrices,  such  as  M,  the  first  element  of  the  eigenvector 
hn  can  be  computed  to  relative  precision,  not  not  just  to  absolute  precision.  Consequently,  if 
Kq  is  computed  as  in  [H],  formula  (196|)  yields  an  with  high  relative  precision. 

5  Numerical  results 

In  this  section  we  present  results  of  several  numerical  experiments.  The  algorithms  for  comput¬ 
ing  the  left  singular  functions  vn  and  singular  values  an  of  Ca (the  operator  defined  in  (l24l)l 
were  implemented  in  FORTRAN  77,  using  double  precision  arithmetic,  and  compiled  using 
GFORTRAN. 

In  Figures  |T]  [2]  and  [3]  we  present  examples  of  left  singular  functions  of  the  operator  Ca,b , 
where  a  =  1  and  b  =  1.1,  b  =  10  and  b  =  100000  respectively. 

In  Figure U] and  Tabled]  we  present  the  singular  values  an  of  the  operator  Ca,b,  for  several 
ratios  7  =  6/a;  an  depends  only  on  7  and  n  (see  Lemma  [Toll .  In  tabled]  we  present  several 
singular  values  smaller  than  10~1000;  the  Fujitsu  compiler  with  quadruple  precision  was  used 
in  this  experiment. 

In  Figure [5] and  Table [3]  we  present  the  eigenvalues  of  the  matrix  M ^  defined  in  (1661). 

In  Figure  |6]  and  Table  [4]  we  present  for  several  combinations  of  7  and  n  the  position  of 
the  last  significant  coefficient  r]n  in  the  expansion  defined  in  (|49f) .  that  is  larger  in  magnitude 
than  e  =  10-16.  In  numerical  computations,  the  vectors  r]even,m  (defined  in  (11001)1  and  ijodd’m 
(defined  in  (11031)1  are  truncated  at  about  half  this  number  (see  Remark  14.  II). 

In  figure  [7]  we  present  the  CPU  time  required  for  the  computation  of  the  expansion  of  the 
101-st  right  singular  function  uioo  of  Tg-y,  for  varying  7;  The  experiment  was  performed  on  a 
ThinkPad  X230  laptop  with  Intel  Core  i7-3520  CPU  and  16GB  RAM. 


Figure  1:  Left  singular  functions  of  £17. i-  Figure  2:  Left  singular  functions  of  £170- 
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Figure  3:  Left  singular  functions  of 
A, looooo- 


Figure  4:  Singular  values  an  of  Ca>b,  with 
7  =  b/a. 


Figure  5:  Magnitude  of  the  eigenvalues  of 
the  matrix  M ^  defined  in  (IHHl). 
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Figure  6:  The  position  of  the  last  signif¬ 
icant  coefficient  larger  in  magnitude  than 

io-16. 
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Figure  7:  CPU  time  required  for  computing  the  expansion  of  the  101-st  left  singular  function 
of  C ii7,  as  a  function  of  7.  The  experiment  was  performed  on  a  ThinkPad  X230  laptop  with 
an  Intel  Core  i7-3520  CPU  and  16GB  RAM. 


Table  1:  Singular  values  an  of  Cayb 


n 

7  =1.0E+01 

7  =1.0E+04 

7  =1.0E+07 

7  =1.0E+10 

0 

1. 0235617  +  00 

1.5568717  +  00 

1.67320 E  +  00 

1.7159517  +  00 

1 

3.0987817  -  01 

1.1228817  +  00 

1.4310717  +  00 

1.56644 E  +  00 

2 

8.3956717  -  02 

7.39927 E  -  01 

1.1487017  +  00 

1.36792 E  +  00 

3 

2.23263 E  -  02 

4.7317317  -  01 

8.9221517  —  01 

1.1606417  +  00 

4 

5.9002017  -  03 

2.9969717  -  01 

6.82645 E  -  01 

9.6834417  -  01 

10 

1.94760 E  -  06 

1.8633617  —  02 

1.28322 E  -  01 

2.9645617  -  01 

20 

3.0080517  -  12 

1.7796717  —  04 

7.7003417  -  03 

3.9511317  -  02 

40 

7.1141517  —  24 

1.60942 E  -  08 

2.7486217  -  05 

6.9538917  -  04 

100 

9.3435917  -  59 

1.1817917  -  20 

1.24105 E  -  12 

3.7635017  -  09 

200 

6.81449 E-  117 

7.0456617  -  41 

7.0878917  -  25 

6.2632517  -  18 

300 

4.1988017  -  61 

4.0463717  -  37 

1.0419017  -  26 

400 

2.50198 E  -  81 

2.3097717  -  49 

1.7330517  -  35 

500 

1.4908117  -  101 

1.31842 E  -  61 

2.8825417  -  44 

600 

8.8829117  -  122 

7.52539 E  -  74 

4.7943717  -  53 

700 

5.2927517  -  142 

4.2953617  -  86 

7.97413 E  -  62 

800 

2.4517017  -  98 

1.32627 E  -  70 

900 

1.39937 E-  110 

2.2058517  -  79 

1000 

7.9872417  -  123 

3.6687817  -  88 
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Table  2:  Examples  of  singular  values  an  smaller  than  10 


-1000 


7 

n 

(%-n 

1.1.E  +  0 

520 

8.70727 'E  -  1002 

1.0E+  1 

1721 

3.6693415  -  1001 

1.0E  +  2 

2797 

5.2996115  -  1001 

1.0S  +  3 

3872 

5.71146 E  -  1001 

1.015  +  4 

4946 

9.4419115  -  1001 

1.0S  +  5 

6021 

8.8974815  -  1001 

Table  3:  Eigenvalues  of  M ^  defined  in  (l66li 


n 

7  =1.0E+01 

7  =1.0E+04 

7  =1.0E+07 

7  =1.0E+10 

0 

—  1.3708115  +  00 

— 7.6814715  +  02 

— 6.8566715  +  05 

-6.58542 E  +  08 

1 

— 4.9931015  +  00 

— 1.2439215  +  03 

— 8.7438615  +  05 

— 7.6083615  +  08 

2 

— 1.2217015  +  01 

— 2.1239415  +  03 

— 1.2050615  +  06 

-9.35829 E  +  08 

3 

— 2.3056115  +  01 

— 3.4392415  +  03 

— 1.6890115  +  06 

— 1.1876915  +  09 

4 

-3.75087£  +  01 

— 5.1952015  +  03 

— 2.3319215  +  06 

— 1.5194715  +  09 

10 

—2.0010215  +  02 

— 2.4969415  +  04 

— 9.5753815  +  06 

— 5.2421315  +  09 

20 

— 7.6014715  +  02 

— 9.3087715  +  04 

— 3.4538415  +  07 

— 1.8075915  +  10 

40 

— 2.9641915  +  03 

— 3.6116715  +  05 

— 1.3278215  +  08 

— 6.8586915  +  10 

100 

— 1.8248015  +  04 

— 2.2201415  +  06 

— 8.1404715  +  08 

— 4.1885315  +  11 

200 

-7.26265 E  +  04 

— 8.8342415  +  06 

— 3.2379315  +  09 

— 1.6650715  +  12 

300 

— 1.9843115  +  07 

-7.27238 E  +  09 

— 3.7393415  +  12 

400 

— 3.5246715  +  07 

— 1.2917415  +  10 

— 6.6416715  +  12 

500 

— 5.5045015  +  07 

— 2.0172915  +  10 

-1.03720£  +  13 

600 

— 7.9238115  +  07 

— 2.9039015  +  10 

— 1.4930515  +  13 

700 

— 1.0782615  +  08 

— 3.9515715  +  10 

— 2.0317015  +  13 

800 

— 5.1602915  +  10 

— 2.6531515  +  13 

900 

— 6.5300715  +  10 

— 3.3574115  +  13 

1000 

— 8.0609015  +  10 

— 4.1444715  +  13 
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Table  4:  The  position  of  the  last  significant  coefficient  larger  in  magnitude  than  10 


-16 


n 

7=  1.0E+01 

7  =  1.0E+04 

7  =  1.0E+07 

7=  1.0E+10 

0 

50 

1502 

43890 

1282730 

1 

51 

1547 

45085 

1318363 

2 

54 

1580 

45902 

1341354 

3 

57 

1611 

46583 

1359605 

4 

58 

1638 

47194 

1375474 

10 

70 

1788 

50346 

1453028 

11 

71 

1813 

50831 

1464669 

20 

88 

2016 

54968 

1562968 

21 

89 

2037 

55411 

1573411 

40 

122 

2438 

63464 

1762452 

41 

123 

2459 

63873 

1772049 

100 

216 

3602 

86776 

2306888 

101 

219 

3621 

87149 

2315585 

200 

366 

5402 

122682 

3143252 

201 

369 

5419 

123031 

3151357 

300 

7124 

156872 

3937590 

301 

7141 

157209 

3945395 

400 

8804 

190098 

4707868 

401 

8821 

190427 

4715479 

500 

10458 

222686 

5462054 

501 

10475 

223009 

5469531 

600 

12092 

254816 

6204538 

601 

12107 

255137 

6211915 

700 

13710 

286598 

6938044 

701 

13727 

286915 

6945339 

800 

318102 

7664390 

801 

318415 

7671623 

900 

349378 

8384872 

901 

349691 

8392049 

1000 

380468 

9100436 

1001 

380777 

9107569 
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6  Conclusions  and  generalizations 


In  this  paper  we  introduced  effective  algorithms  for  the  evaluation  of  the  left  singular  functions 
and  singular  values  of  the  Truncated  Laplace  Transform  Cafi-  Together  with  the  algorithms 
introduced  in  [12]  for  the  computation  of  the  right  singular  functions,  these  algorithms  conclude 
the  construction  of  the  SVD  of  the  operators  Ca,b- 

As  is  evident  from  Remark  12.31  and  the  more  detailed  discussion  in  [20] .  the  left  singular 
functions  of  Ca,b  are  an  efficient  basis  for  representing  combinations  of  decaying  exponentials 
whose  decay  constants  are  in  the  interval  [a,b\. 
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